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Abstract 

The generation of primordial magnetic seed fields during inflation is studied in a theory derived 
from the one-loop vacuum polarization effective action of the photon in a curved background. This 
includes terms which couple the curvature to the Maxwell tensor. The resulting magnetic field 
strength is estimated in a model where the inflationary phase is directly matched to the standard 
radiation dominated era. The allowed parameter region is analyzed and compared with the bounds 
necessary to seed the galactic magnetic field. It is found that magnetic fields of cosmologically 
interesting field strengths can be generated. 



1 Introduction 

Magnetic fields are found to be associated with nearly all structures in the universe. They are 
observed on stellar upto possible supercluster scales pQ. In some processes such as star formation 
they can play an important role. Equally the physics of cosmic rays indicate the existence of a large 
scale galactic magnetic field. Furthermore, ultra high energy cosmic rays which are most likely of 
extra galactic origin could be used to study the properties of galactic and extra galactic magnetic 
fields |2]. 

As to the origin of the observed magnetic fields there seem to be generally speaking two broad 
classes of mechanisms [lj. On the one hand there are battery-type mechanisms which work on the 
basis of charge separation which leads to a current and finally induces a magnetic field. On the 
other hand there are dynamo type mechanisms which amplify an initial seed magnetic field. The 
former have a coherence scale of the order of the domain associated with the battery-mechanism, 
which in a cosmological context is always smaller than the horizon at the epoch of creation. The 
latter involves magnetic fields whose correlation length is of the order of the region of interest, which 
in the case of a galactic magnetic field, would be a proto galactic scale of the order of 1 Mpc today 
which requires a mechanism to create magnetic fields with rather large coherence scales. A natural 
mechanism for this is provided by the amplification of perturbations of the electromagnetic field 
during inflation. Since quantum perturbations are stretched beyond the horizon during inflation 
and upon leaving the causal domain becoming classical, there is no problem with the coherence 
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length. However, in general, in a background geometry with flat spatial section the resulting field 
strength of the primordial magnetic field after inflation in standard electrodynamics is far too 
small to, for example, serve as a seed field for a potential galactic dynamo explaining the observed 
galactic magnetic field of the order of 10 -6 G. Starting with [3] this motivated an intensive study 
of alternative models of some kind of coupling of linear electrodynamics to either other fields in 
the theory, such as scalar fields [3j or gravity [31 [5] including extra dimensions [6] . Other models 
break explicitly Lorentz invariance considering a non zero photon mass [7] . Recently there has 
also been interest in magnetic field generation within models of nonlinear electrodynamics which 
naturally occur when quantum corrections and self couplings of the electromagnetic field are taken 
into account [8]. Furthermore, there are models within the standard model or its super symmetric 
extensions [9]. In the case of models with curved spatial sections it was shown that even in the 
minimally coupled model of linear electrodynamics strong enough magnetic fields can be created 

Here we are returning to a model where the electromagnetic field is gravitationally coupled 
which first has been proposed in the context of the generation of primordial magnetic fields during 
inflation in [3]. In particular the model under consideration derives form the one loop effective 
action of vacuum polarization in QED in a gravitational background [11]. As shown in [11] this 
leads to birefringence of the electromagnetic wave where the photons have velocities depending on 
their polarization which can exceed the speed of light. This last observation leads to an interesting 
causal structure of these space-times [12] . 



2 Estimating the magnetic field strength 

In [11] it was shown that to lowest order the propagation of a photon in a gravitational background 
is described by the Lagrangian [ITJ [3] 



1 1 

C = —-F F^ v - 

4 ^ Ami 



bRF^F^ + cR^F^F\ + dR^ XK F^F XK + f(V^)(V a F a u ) (2.1 



The expansion parameter is bascially the square of the Compton wave length of the electron which 
enters due to the fact that vacuum polarization effectively gives the photon a non zero "size" due 
to electron positron pair creation. The parameters b, c, d, and / are are free parameters, which, 
however, have been calculated in [11] in the weak gravitational field limit. It was found that the 
coefficient of the last term is of the order e 2 and can therefore be neglected. Moreover, unless / is 
to be chosen much larger than the other coefficients it is expected to be negligible in general, as 
can be appreciated from the equations of motion, 



mi 



bRF^ + C - [R\F Xu - R\F Xll ) + dR XK ^F Xh 



+ (VaV^V^ + R a u V^F Pa ) = 0. (2.2) 

Thus the new type of term that the last term contributes involves higher derivatives which are 
suppressed by a factor ~ 1/<t 2 , where a is a typical scale. Therefore, in the following this term will 
neglected and / = 0. 

The Maxwell tensor is expressed in terms of the gauge potential A^, that is F^ = d^A u — d u A^ 
and the Coulomb gauge is used Aq = 0, diAi = 0. Furthermore, the background cosmology is 
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described by the line element, 

ds 2 = a 2 (rj) ( - drj 2 + dx 2 + dy 2 + dz 2 



(2.3) 



where a(rj) is the scale factor describing a model, in which the inflationary stage is directly matched 
at r] = 771 to the standard radiation dominated era, 



a(rj) 



ai 



r\-2r\\ 

-m 



f] < rji 



v > m- 



(2.4) 



In the following a\ = 1. For (3 = — 1 de Sitter inflation is realized and for —00 < (3 < — 1 the model 
has a stage of power law inflation for r] < rj\. The exponent (3 is related to the equation of state of 
matter, defined by, p = 7/5, by (3 = 3 2 +1 . Note that 7 = takes values between -1 and — |. The 
Fourier expansion of the gauge potential is given by, 



Aj(r],x) 



d 3 k 



£ (A) 

kj 



,(A)t 



ik-x 



where the sum is over the two polarization states and e- are the polarization vectors satisfying, 

kj 

ei • k = 0. Furthermore, the amplitude A k satisfies the same mode equation for both polarization 
states. Thus the index (A) is suppressed in Ak- 
in Fourier space in the background model f|2.3[) equation (|2.2p results in 



(2.5) 



Fifa)A fc + F 2 (rj)A k + F 3 ( V )k 2 A k = 0, 



(2.6) 



where a dot indicates and 
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/*i 



-2(/3+l) 



-2/3-3 
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7/ 2 m 2 
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/xi = j9 66(/9 - 1) + c(/3 - 2) - 2d 
H2 = -2(/3 + l)/ii 

^3 = /j[6&(/? - 1) + c(2^ - 1) + 2d/3 



(2.7) 



The standard quantization procedure requires the corresponding action of the field to be diagonal. 
This can be achieved by using the canonical field ^ defined by 



* = F 1 *A k 



(2. 



which satisfies the mode equation, 



q>" + pq; = 0, 



(2.9) 
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where a new dimensionless variable z = —krj has been defined and ' = 4-. Moreover, 



P = - 



-4/3-6 



1 + KoZ-W+V 



1 

+ -x 



K 3 Z 



-2/3-4 



1 + KiZ -2{P+l) 



+ 



(2.10) 



and 



where kq 



2 2 
f i 2 K 



m, 



Hi 



-2 



K 2 = HlKQ 
k \ 208+1) 



k 3 = (2/3 + 3)/i 2 K 



K 4 = /X 3 K 



(2.11) 



In deriving these expressions the maximally amplified (comoving) wavenumber k% has been defined 
by k\ = T—y H\ is the value of the Hubble paramter at the beginning of the radiation dominated 
stage at r}\. It is related to k\ by k\ ~ H\. The spectrum of the resulting magnetic field is 
determined by calculating the Bogoliubov coefficients which connect the "in" and "out" vacua. 
The relevant solutions of the mode equation (|2.9|) are those on superhorizon scales, z <C 1, in which 
case equation (|2.9[) reduces to 



V' + {tiz- 2 + 6)^ = 0, 



(2.12) 



where 
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66Q9 - 1) + c(2/3 - 1) + 2d(3 
6b(/3 - 1) + c((3 - 2) - 2d ' 



(2.13) 



For /3 = — 1, which describes de Sitter inflation, £i = and £2 = 1 the solution is a plane wave 
which was also noted in [HJE]. Furthermore, j3 = — 2 also gives £1 = 0, but £2 = i86+4c+2rf • Solving 
equation (|2.12p during power law inflation, j3 < — 1 and [3 7^ —2, results in the following solution 

(2) 

in terms of the Hankel function of the second kind, H{, (x), 



* (I) = V21^ 2)( ^ Z) 



where 



+ 



(2.14) 



which gives the correctly normalized incoming wave function for 7/ — > — 00 for £2 > . During 
the radiation dominated stage the additional curvature terms can be neglected and thus the mode 
equation (|2.12p simplifies to that of a free harmonic oscillator which is solved by the superposition 
of plane waves, 



(RD) 



1 



(2.15) 



where zi = fc|r/i| and c± are the Bogoliubov coefficients. Since the aim here is to calculate the 
magnetic field energy spectrum at galactic scale which re-enters during the radiation era it is enough 
to only consider the radiation dominated stage. Matching the solutions of the gauge potential and 
its first derivatives at rj = rji on superhorizon scales determines c+ and c_. In particular, using the 
small argument limit of the Hankel functions [13J |c_| 2 is found to be, for (3 7^ — | 



I c 1 2 ~ 
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(2.16) 
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where it was used that in the approximation used here, Fi(j]i) ~ \i\ (jjpj^j ■ In the case (3 = — § 
the limiting behavior of the mode function on superhorizon scales leads to a a divergent factor 
In 2 ^y/^2j^j in |c_| 2 . Therfore, this case will not be considered further. Including both polarization 
states the total spectral energy density of the photons is given by (cf., e.g., [H] ) 

, \ d P /A;\ 4 |c_| 2 . . 

a log k \a J 7r z 

Since the electric field decays rapidly due to the high conductivity of the radiation dominated 
universe, the spectral energy density (|2.17|) gives a measure of the magnetic field energy density, 
ps, which results in the commonly used ratio of magnetic field over background radiation energy 
density r = ^ = ^ [3], where the density parameter of radiation is given by 7 = (-^) (^-) 4 , 
for/3^-2,-1,-1, 7 



2[I»] 2 (I \ 2 / mA 2 /6 VVA N " 2 ' 



V2--j U?J If J U' (2 ' 18) 



where Mp is the Planck mass. An expression similar to (|2.18p is obtained when the magnetic field 
energy density is calculated using the two point function of the magnetic field, (Bi(k)B*(k')). The 
form of the magnetic field spectrum (|2.18p imposes the constraint v < | which implies the range 
for (3 given by —3 < (3 < — 1 taking into account the constraint from power law inflation. 

An initial magnetic field with strength B s ~ 10~ 20 G could seed the galactic magnetic field 
assuming that in addition there is a galactic dynamo operating [15J. This leads to r ~ 10" 37 at 
a galactic scale of order 1 Mpc. The former estimate does not take into account the presence of 
a cosmological constant which reduces the minimal magnetic field strength to B s ~ 10 -30 G and 
r ~ 10~ 57 [16] . In order to directly seed the magnetic field, r has to be of the order of 10~ 8 . In the 
following r (cf. equation ([2.18]) ) is calculated at the galactic scale wq = 10 _14 Hz, corresponding 
to a length scale of 1 Mpc. Furthermore the physical frequency corresponding to the maximally 

amplified wave number, k\, today, is given by lui(t]q) ~ 6 x 10 11 (^jj^J^ Hz [6j- The expression 

for r (cf. equation (|2.18j) ) depends on the parameters b, c and d in addition to the parameters 
characterizing inflation (3 and the Hubble parameter at the beginning of the radiation dominated 

stage -p^-. The corresponding reheat temperature T% is given by ~ J An upper limit on r 

can be derived by using the constraint 

Mi fe) > 1, (2-19) 

consistent with the approximation used to derive equation (|2.12p . This yields to, for (3 ^ — 2, — |,— 1, 
at oj g = 10" 14 Hz 

W(^) = 1(T 79 + 52 ' [I»] 2 Q - „) 2 (|) U . (2.20) 

For simplicity assuming that the parameters determining the curvature terms are all of the same 
order, so that b ~ c ~ d, the maximal possible ratio r max is independent of these parameters, since 
£2 in this case is given by £2 = 7/3-10 • fig ure [H (left) t ne contour lines for log 10 r max are shown 
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b=c=d 



b=c=d 




Figure 1: Left: The contour lines for the maximum value of the logarithm of the ratio of magnetic 
to background radiation energy density is shown, that is log 10 r max , for b = c = d. The values of 
log (^-jq^J correspond to reheat temperatures between 0.1 GeV and 10 19 GeV. The numbers within 
the graph refer to the value of log 10 r max along the closest contour line. 

Right: The contour lines for the minimum value of the logarithm of the parameter 6, that is 
lo 

§io ^minj are shown in the model b — c — d. The numbers within the graph refer to the value of 

lo 

gio ^min along the closest contour line. 



in this case as a function of (3 and log 10 (^jj^J- The constraint on fj,\ (cf. equation (12,19h ) implies 
in the case b ~ c ~ d a constraint on the parameter b, 



b > K 



10 



-45 



0(7/3-10) \M P J 



(2.21) 



log 6min is shown in figure Q] (right). 

In [IT] a strong bound on primordial magnetic fields created before nucleosynthesis has been 
derived. The bound is due to the conversion of magnetic field energy into gravitational wave energy. 
Translating the bound on the magnetic field strength given in [17] into a bound on the magnetic 
field to background radiation energy density ratio r results in 



r GW * 2 x 10- 61+52 ^ GW ^ (H 
at the galactic scale used here, A = 1 Mpc, and 

Mgw 

Therefore, the requirement r max < tqw implies an upper limit on , namely 



2" U T | - 



(2.22) 



(2.23) 



log 



10 



= 9 + 2 lo §io 



2l-^r(|-^) (i 2 y 



(2.24) 
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so that the allowed range is given by yj]~J — yjlpj ■ This is plotted in figure [2] for small values 

of (3. As can be appreciated from figure [2] the maximally allowed value for log 10 (l7p) i s always 
much bigger than 1 for (3 < — 2 which is the relevant range here. This can always be satisfied since 
it is assumed that H% < Mp. 



b=c=d 




Figure 2: The maximal value log 10 ( j^- J (cf. equation (I2.24|) ) allowed by the bound from 

V p / max 

gravitational wave production is shown for the model with b = c = d for h$ = 0.73. 



3 Conclusions 

Primordial magnetic field generation has been investigated in a model where the electromagnetic 
field is coupled to various curvature terms which is motivated by the form of the one loop effective 
action of vacuum polarization in QED in a gravitational background [3]. Here the resulting 
magnetic field spectrum has been calculated explicitly by matching a stage of power law inflation 
directly to the standard radiation dominated phase, solving the mode equation and finding the 
appropriate Bogoliubov coefficient. The ratio of magnetic field energy density over background 
radiation energy density r calculated at a galactic scale of 1 Mpc has been employed to compare 
the model with the minimal required values necessary to seed the galactic dynamic field either 
directly or relying on the mechanism of a galactic dynamo. For simplicity it was assumed that the 
coefficients of all additional curvature terms are of the same order. A constraint on the parameters, 
which is part of the approximation used to derive the mode equation on superhorizon scale, leads 
to a maximum value of the magnetic field strength or equivalently the ratio r which is a function 
of (3 and Here j3 is the exponent characterizing inflation and H\ is the value of the Hubble 
parameter at the beginning of the radiation dominated era. As can be seen from figure [T] there is a 
region in parameter space where the resulting magnetic field strength is larger than B s > 10 -20 G 
corresponding to r > 10~ 37 , namely, (3 < —2.4, -|p > 1CP 18 bounded by the corresponding contour 
line. De Sitter inflation corresponds to (3 = — 1. In this case there is no significant magnetic field 
generation since the mode functions during inflation as well as during the radiation dominated 
stage are plane waves. 

Finally, it has been checked that the resulting maximum magnetic field strength satisfies the 
bound due to gravitational wave production [17] . Imposing this bound leads to a maximal value of 

which, however, is always much larger than the maximal value considered here, as can be seen 
from figure [2J In summary, there is a range of parameters (3 and for which magnetic fields are 
generated that are strong enough to explain the galactic magnetic field. 
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